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[57] 



ABSTRACT 



An improved Gaussian minimum shift keying (GMSK) 
carrier tracking loop operating at baseband takes advantage 
of the orthogonality of precoded data GMSK signals and 
Laurent filtering to provide a carrier phase error signal 
generated at baseband for carrier phase derotation of the 
received GMSK signal. The carrier tracking loop also pro- 
vides demodulated data estimates with performance equal to 
that of a GMSK serial demodulator. The tracking loop uses 
data directed feedback to improve noise rejection, but still 
has fast acquisition by operating at baseband. 

13 Claims, 3 Drawing Sheets 




62 



REAL 
INPHASE 
h Q (-t) 



^•60 



COMPLEX TO 
REAL AND 
IMAGINARY 
CONVERTER 



64 



I 



IMAGINARY 
f\ QUADRATURE 

h 0 (-t) 



58 



52 



L e -i<o c t 



NCO > 
90 



66 



AI 



t = 2kT 



t = {2k+irrY 









I 



70 



A 



68 



Obi 

t = 2kT \ 



JL 



76 



A 



72 



t=(2k*irr\^ e BQ v O / J ( 




Q 

80 



0© 





r 84 








even 



PRECODED GMSK CARRIER TRACKING LOOP 



10/13/2004, EAST Version: 1.4.1 



U.S. Patent 



Nov. 14, 2000 



Sheet 1 of 3 



6,148,040 



NF 
DA 
SOU 


IZ 

TA 

RCE 




d(t) 



18 



I 

22 



▼ r 



20 



DATA 
PRECODER 



a(t) 



10 



I 



INTEGRATOR 
^26 

g(t) 



GAUSSIAN 
FILTER BT b 



TRANSMITTER iMpDULAJOR 



G(t) 



28 



PHASE AND 

CARRIER 
MODULATOR 
I 



h 

24 



CARRIER 
REFERENCE 



•30 




s(t) 



II 
II 
II 
II 



GMSK TIMING RECOVERY LOOP 



II 

II 
II46 

l^ 



DIGITAL 
TRACKING 
TRANSITION 
LOOP 



To 



c H (t) 



HARD 
LI MITER 



48 



n j 



r 



42 



T(t) 



■50 



DATA 
SINK 



d(t) 



DATA 
DETECTOR 

C 



R 0 (t) 



40 



DEMODULATOR 



38 



GMSK TIMING 
ACQUISITION LOOP 



36 



PRECODED 
GMSK 
CARRIER 
TRACKING 
LOOP 



R(t) 



8 f 



44 



32 



34 



GMSK CARRIER 
PHASE ACQUISITION 
LOOP 



GMSK CARRIER 
A FREQUENCY ACQUISITION 
LOOP 



RECEIVER 



r. 



12 



PRECODED GMSK COMMUNICATION SYSTEM 

FIG. 1 



10/13/2004, EAST Version: 1.4.1 



U.S. Patent 



Nov. 14, 2000 Sheet 2 of 3 




66 





t = 


r fSZ 




REAL 






INPHASE 




hg(-t) 


t = 


(2k+1)T 



r 



60 



I 



' e AQ 









I 



70 



A 

A 



n 



•68 



COMPLEX TO 
REAL AND 
IMAGINARY 
CONVERTER 



R c (t) 



64 



I 



IMAGINARY 
f\ QUADRATURE 
h 0 (-t) 



t = 2kT \ L— 
W2 



58 



76 



A 



.-(j9 e +e c) 



r 



52 



e -JW c t 



NCO "N 
90 



t = (2K+irr\ e BQ vV 



LOOP FILTER 



88 



©0 




PRECODED GMSK CARRIER TRACKING LOOP 

FIG. 2 



10/13/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 14,2000 Sheet 3 of 3 6,148,040 



13 

12 

BT b = 0.2 

11 
10 
9 

CARRIER 
TRACKING 

PHASE 0 

JITTER 
(DEGREES) 7 

6 

5 

4 

3 

2.90 
IRRIDUCIBLE 
TRACKING 2 f 
JITTER 












































\ 
















































































































\ 
























































































































































































































































) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 



BIT SNR dB 
CARRIER PHASE JITTER PLOT 



FIG. 3 



10/13/2004, EAST version: 1.4.1 



6,1- 

1 

PRECODED GAUSSIAN MINIMUM SHIFT 
KEYING CARRIER TRACKING LOOP 

REFERENCE TO RELATED APPLICATIONS 

The present application is related to assignee's copending 
application entitled "Gaussian Minimum Shift Keying 
(GMSK) Precoding Communication Method" Ser. No.: 
09/390,966, filed Sep. 07, 1999, by the inventors G. Lui and 
K Tsai, and related to assignee's copending application 
entitled "Digital Timing Recovery Loop for GMSK 
Demodulators" Ser. No. 09/307,231, filed May 07, 1999, by 
the inventors T. Nguyen, J. Holmes, and S. Raghavan, both 
of which are here incorporated by reference as there fully set 
forth. 

STATEMENT OF GOVERNMENT INTEREST 

The invention was made with Government support under 
contract No. F04701-93-C-0094 by the Department of the 
Air Force. The Government has certain rights in the inven- 
tion. 

FIELD OF THE INVENTION 

The invention relates to the field of continuous phase 
modulation communication systems. More particularly, the 
present invention is related to Gaussian minimum shift 
keying carrier tracking loops for use in combination with 
data preceding. 

BACKGROUND OF THE INVENTION 

Communication systems have long transmitted digital 
signals using various carrier modulation techniques. The 
spectrum of a digital signal can be controlled and made 
compact by envelope filtering or phase domain filtering. An 
efficient phase domain filtering approach controls the signal 
spectrum by frequency modulating the filtered signal onto a 
carrier frequency to form a continuous phase modulated 
(CPM) signal. Because the CPM signal has a constant 
envelope, a power amplifier can be operated at maximum 
output power without affecting the spectrum of the filtered 
signal. Gaussian minimum shift keying (GMSK) is a form of 
continuous phase modulation. GMSK uses CPM signals 
with a constant signal envelope and a spectrum that can be 
made compact with the appropriate choice of the signal 
bandwidth bit time product (BT) product. 

An M-ary GMSK signal is defined by its complex enve- 
lope described in terms of symbol energy E, symbol period 
T, carrier phase 0 C and phase modulation 6(t) using a 
modulation index h. Input data is formatted into data sym- 
bols prior to carrier modulation and transmission. The data 
formatting may be non-return to zero (NRZ) formatting. 
Equally probable NRZ data symbols belong to an M-ary 
alphabet of symbols having the symbol time T. The M-ary 
symbols are used to phase modulate a carrier reference. The 
GMSK phase response 8(t) originates from a Gaussian filter 
response g(t) of a Gaussian smoothing filter with a single 
sided 3 dB bandwidth B, truncated to an intersymbol interval 
duration L, that is a memory truncation length L. The GMSK 
Gaussian filter with a memory truncation length L of a 
GMSK signal is defined by the BT bandwidth bit time 
product, where B is the single sided 3 dB filter bandwidth in 
hertz. The Gaussian filter with a small BT product, has a 
memory length L equal to 1/BT The Gaussian filter response 
g(l) used to phase modulate the carrier by a phase modulator 
having a modulation index h. In general, lowering the 
modulation index h while keeping the BT product constant 
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will further reduce the spectral occupancy. The intersymbol 
memory length L is the number of elapsed symbol periods 
for the GMSK signal to accrue a complete phase change 
amount due to a single input symbol and hence represents 

5 the memory of the GMSK signal. The phase modulated 
GMSK signal is transmitted to GMSK receiver for commu- 
nicating the input data stream. 

The GMSK receiver demodulates the received GMSK 
signal into a demodulated signal that is in turn passed 

10 through Laurent filters providing filter signals fed into a 
Viterbi decoder for providing a estimate of the input data 
stream. 'The Laurent filters are applied to an accumulated 
phase at a current bit time. The current bit has signal 
components extending over Lbit periods. A typical coherent 

15 receiver for M-ary GMSK signal is based on a pulse 
amplitude modulation (PAM) representation of CPM signals 
using Laurent filtering, and employs the Viterbi algorithm to 
optimally demodulate symbol sequences. In demodulating 
M-ary GMSK signals using the Viterbi algorithm, a differ- 

20 ential decoder has been necessary to resolve data bit ambi- 
guities while providing a nominal bit error rate (BER) that 
is desirably as small as practicable. A Viterbi algorithm 
typically employs a sliding window in the demodulation 
process where the width of the sliding window represents 

25 the demodulation memory or delay. The surviving state 
sequence produced by the sliding window Viterbi algorithm 
at stage o depends on all the demodulated symbols d„(l). The 
intrinsic data dependency of the survivor sequences a„(t) 
disadvantageous^ requires a differential decoder operation 

30 in the receiver when deciding on the actual demodulated 
symbol from successive survivors of the Viterbi algorithm 
resulting in a differential bit error rate degradation. 

In the related application, Nguyen et. al., an improved 
GMSK timing recover loop offers closed loop generation of 

35 a data timing signal at a baseband frequency. The improved 
GMSK timing recovery loop enables recovery of the trans- 
mitted data using the baseband data timing signal x(t) with 
high accuracy at low bit signal to noise ratio (BSNR) and at 
a small BT product, and has the advantage of negligible loss 

40 due to non-random data patterns. Another advantage asso- 
ciated with GMSK timing recovery loop is the adoption of 
the well known digital transition tracking loop (DTTL) used 
in M-ary PSK systems with a modification of adding a hard 
limiter in closed loop control of the data timing signal. 

45 The GMSK system includes the transmitter modulator 
and the receiver demodulator between which is transmitted 
the GMSK signal. The demodulator includes a carrier track- 
ing loop for providing a GMSK demodulated received signal 
R 0 (t) and the bit timing recovery loop for providing the bit 

50 timing signal x(t). The carrier tracking loop preferably 
employs reverse modulation. The GMSK timing recovery 
loop performance employs the hard limiter adjusted by a bit 
timing error signal for improved insensitivity to the values 
of BT while operating at low BSNR. The GMSK timing 

55 recovery loop takes advantage of the observation that the 
cosine of the baseband GMSK signal has zero crossings at 
multiples of the bit duration. The hard limiter is used to 
create the clocking signal for the NRZ data stream that has 
the zero crossings at multiples of the bit duration. The digital 

60 transition tracking loop is then used to track the zero 
crossings of the NRZ data stream clocking signal from the 
received demodulated GMSK signal, and the bit timing 
signal x(t) is then generated by the DTTL with less jitter for 
improved data detection. In the GMSK timing recovery 

65 loop, the hard limiter is adjusted by the bit timing error 
signal x e (t) to reduced jitter in clock sampling of the NRZ 
data stream. Hence, the digital transition tracking loop tracks 
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the adjusted zero-crossings of the NRZ data stream, and the Another object of the invention is to provide carrier phase 

reduced jitter bit timing signal x(t) is then generated for tracking operating baseband frequencies of a demodulated 

accurate data sampling and detection. Significantly, this received signal. 

timing recovery loop is operated at baseband and is a Yet aQOther ob - c( of ^ mvention is t0 provide c^er 
preferred improvement to the GMSK receiver. 5 faase ^ 0 atin al baseb and frequencies to con- 
In the related application, Lm et. al., a data preceding ^ r consun3ption . 
algorithm is implemented prior to modulation in the trans- - — .— - ; — . — ~~ " — ; -_ _ WOT - 
miner to substantially improve the resulting BER perfor- (The present mvention is duected to an improved GMSK ^ 
mance of the continuous phase modulated (CPM) transmit- carncr ***ing 1°°P operatmg at baseband while taking 
ters and receivers, such as the Gaussian minimum shift 1Q advantage of the orthogonality of precoded data GMSK 
keying (GMSK) transmitters and receivers without the use signals and the absolute phase of the PAM filtered signals, 
of differential decoders while preserving the spectral occu- As a precondition, the transmitted data sequence d(t), is 
pancy the GMSK signals. The preceding algorithm encodes precoded in the transmitter and communicated by Gaussian 
the source NRZ data symbols prior to the GMSK modula- filtering and phase modulation of a carrier reference f c lto 
tion so that the cumulative phase of the precoded symbols provide a transmitted GMSK signal to a receiver. Trie 
becomes the absolute phase of the data symbols in the signal 35 GMSK carrier tracking loop carrier and phase demodulates 
phase trellis of the Viterbi algorithm. The preceding algo- the received GMSK signal R(t) into a baseband demodulated 
rithm offers performance improvement for M-ary coherently received~R^(t) signal that is then separa ted in to real inphase 
demodulated GMSK signals. C(l) and imaginary ^quadrature "phase (Q) components both 
Precoding improves the BER performance for coherent subjected to Laurent filtering to provide rough estimates of 
demodulation of the M-ary GMSK signals implemented 20 the baseband data signal 0 f encoded data. Orthogonal 
using a pulse amplitude modulated signal subject to the switching and hard limiting is used to generate data A„ and 
Viterbi algorithm. The preceding algorithms encodes the g est i ma tes, and of respective Q and I components at 
source NRZ data symbols d„(t) prior to the GMSK modu- baseband . ^ ^ and fi„ data estimates and 6, and 9 G phase 
lation so that the cumulative phase of the precoded symbols em)r ^ fc cress-multiplied together for generating a 
d„(t) is identical to the exact phase of the source NRZ 25 u • . A , ■ u . i • ♦ u u j 
symbols at every stage of the Viterbi algorithm. In the !£ ase erro [ m ' lose < loo P carner trackm S at baseband ' 
Viterbi algorithm, the preceding process produces a set of P 1 ^ u and f ° her ^vantages will become more apparent 
survivor sequences for estimating the original data bit with- fr °^ [ he followin S detailed description of the preferred 
out the use of differential decoding. The Gaussian filter can embodiment. 

be expressed mathematically, and the Laurent mathematical 30 BRIEF DESCRIPTION OF THE DRAWINGS 
expansion dictates the matched filter bank. Without 

preceding, the Gaussian filter creates phase ambiguities that FIG. 1 is a block diagram of a precoded GMSK commu- 

are resolved by differential decoding. Because the precoded nication system. 

symbols have the same statistics as the source symbols, the FIG. 2 is a block diagram of a precoded GMSK carrier 

transmit spectrum of the GMSK signal is preserved while tracking Loop 

eliminating the need for differential decoding. Depending 35 „ T „ - . . - . , . , .. . , 

upon the channel bit error rate in operation, The precoding ™" ? * P ^ " g P ^ m ^ 

method will render a signal to noise ratio (SNR) improve- as a ^ nch0D of blt signal-to-noise ratio, 

ment of 3 dB over the same modem that demodulates GMSK DETAILED DESCRIPTION OF THE 

signabwuhom preceding PREFERRED EMBODIMENT 

The carrier phase demodulation of the received signal has 4U 

long been performed at high intermediate frequencies (IF) An embodiment of the invention is described with refer- 

causing squaring signal losses and consuming high power ence to the figures using reference designations as shown in 

when demodulating at the IF frequencies. The squaring the figures. Referring to FIG. 1, a precede GMSK commu- 

losses disadvantageous^ increase the bit error rate. The nication system includes a GMSK transmitter 10 and 

techniques used for carrier phase synchronization usually 45 receiver 12 for communicating a GMSK signal S(t) from a 

require squaring or Costas loops with losses due to squaring, transmitting antenna 14 to a receiver 16 antenna providing 

and self noise due to intersymbol interference (ISI) with an intermediate frequency (IF) received signal R(t) received 

high BER. These carrier tracking loops do not perform well by a GMSK receiver 12. The GMSK signal S(t) is subjected 

in the presence of non-random data patterns where the to noise and interference, both not shown, during commu- 

discrete components for carrier recovery may vanish. 50 nication of the signal. The transmitter 10 includes a data 

A reverse modulation method may be used in carrier source 18 that may be a non-return to zero (NRZ) formatted 

phase tracking loops operating at high intermediate frequen- data source providing a data stream d(t) to a data precoder 

cies. The reverse modulation method works very well with 20 providing in turn a sequence of precoded symbols a(t) 

differentially encoded data. However, when used with the where *(0 fc defined as a^O-d^Od^-l)" having bit 

precoded data, the tracking performance becomes unstable 55 duration of T for the nth bit with a bit rate of R fc . The 

and sensitive to the loop gain. Additionally, when operating precoded symbols a(t) are communicated to a modulator 22 

at high intermediate frequencies, more power is disadvan- tnat includes a Gaussian filter 24 having a predetermined 

tageously consumed. While prior GMSK systems have used bandwidth bit-time product BT fc . The Gaussian filter 24 

preceding to avoid receiver differential decoding, the pre- receives the encoded symbols a(t) as volt pulses to 

coded data absolute phase characteristics have not been used 60 generate respective Gaussian filter pulse responses g(t) 

for baseband operation of a carrier tracking loop. These and overlapped by memory length L and communicated to an 

other disadvantages are solved or reduced using the inven- integrator 26 providing an accumulative Gaussian filter 

tion. response signal G(t) that is in turn communicated to a phase 



modulator 28 that phase modulates a f c carrier reference 30 
65 by the accumulative Gaussian filler response signal G(l) to 
An object of the invention is to provide carrier phase provide the GMSK signal S(t). Gaussian filter responses g(t) 



SUMMARY OF THE INVENTION 
>ct of the invention is to provide car 
tracking using data precoded GMSK signals. for the symbols a(t) are superimposed over the time period 
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of the memory truncation length L and communicated to the 
integrator 26 when generating a complete Gaussian filter 
response G(t) communicated to the phase modulator 28 
providing the GMSK signal S(t) as a continuous phase 
modulated (CPM) signal. The phase modulator 28 is defined 
by a modulation index h. The respective Gaussian pulse 
responses g(t) of the Gaussian filter 24 is a function of the 
BT product and the truncation memory length L. The phase 
modulation <(>(t) is equal to the jihG(t). The accumulative 
Gaussian filtered response G(t) is phase modulated by the 
phase modulator 28 to generate the CPM GMSK signal S(t) 
having a constant envelope. The receiver 12 receives the 
transmitted GMSK signal S(t) as the received signal R c (t) is 
a function of the transmitted signal and a noise component. 



signal R 0 (t) is a function of the received signal R(t) that can 
expressed as a complex baseband signal. 



15 



The function of the receiver 10 is to generate an estimate 
d(t) of the original data stream d(t). The receiver 12 includ- 
ing a demodulator 32 receives the received signal R(t) for 
generating a baseband signal R G (t) for data recovery. The 
received signal R(t) is also communicated to a GMSK 
carrier phase acquisition loop 34 providing a carrier phase 
estimate 0 C) to a GMSK carrier frequency acquisition pro- 
viding a carrier frequency estimate f c , and to a GMSK 
timing acquisition loop 38 providing an initial tracking 
timing signal z 0 . A data detector 40 generates the estimated 
data d(t) communicated to a data sink 42. The data detector 
40 receives the received baseband signal R 0 (t) from a 
precoded GMSK carrier tracking loop 44 and a data timing 
x(t) signal from a GMSK timing recovery loop 46 preferably 
including a digital tracking transition loop 48 and a hard 
limiter 50. The digital tracking transition loop 48 provides 
the data timing signal T.(t) to the data detector 40 for 
recovery data estimation. The data timing signal x(t) is 
received by the hard limiter 50 providing a hard clocking 
signal Cyy(t) to the digital tracking transition loop 48 in 
closed loop control for rapidly stabilizing the data signal x(t) 
at baseband for clocking the data detector 40. The data 
detector 40 is preferably a trellis receiver having Viterbi 
decoding and Laurent matched filtering. The data detector 
40 preferably includes a Laurent filtering, data sampling and 
Viterbi decoding for providing the data estimate d(t). The 
received baseband signal R 0 (t) is a representation of the 
output of the integrator 26. The timing recovery loop 46 in 
combination with the data detector 40 is an inverse function 
of the data precoder 20, Gaussian filter 24 and integrator 26, 
to provide the estimated data d(t). Loops 34, 36 and 38 may 
be of conventional designs. The invention is directed to the 
precoded GMSK carrier tracking loop 44 shown in detail in 
FIG. 2. 

Referring to FIGS. 1 and 2, and more particularly to FIG. 
2, the improved precoded GMSK carrier tracking loop 44 
receives the carrier frequency estimate f,. and carrier phase 
estimate Q c to provide a carrier demodulated and phase 
adjusted baseband received signal R 0 (t). The carrier signal f c 
is received by a carrier receiver 52 to provide a carrier 
demodulation signal to a carrier demodulator 54 that 
demodulates the received signal R(t) into a carrier demodu- 
lated baseband received signal R c (t) that is then communi- 
cated to a phase demodulator 56 that removes the carrier 
phase 0 C and a phase error 0,. from a phase receiver 58 to 
then provide the carrier and phase demodulated baseband 
received signal R 0 (t). The demodulated baseband received 
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The demodulated baseband received signal R 0 (t) is 
received by a complex to real and imaginary converter 60 to 
provide a real inphase Q component and an imaginary 
quadrature I component respectively communicated to a 
fundamental PAM real inphase filter h 0 (t) 62 and to a 
fundamental PAM imaginary quadrature filter b<,(t) 64 to 
respectively provide a 0, PAM filtered signal and a 0 Q PAM 
filtered signal. The 0 7 and Q Q PAM filtered signals from the 
PAM filters 62 and 64 are baseband signals having time 
varying positive and negative phase values corresponding to 
positive and negative pulses of the input data precoded 
sequence a(t). 



25 



a 0 „ = J{n)d nt J(n) 



evenl 
odd J 



The 0, PAM filter signal from the PAM real filter 62 is 
communicated to a 2kT sample and hold switch 66 provid- 
ing an even sampled phase 6 A/ and communicated to a 
(2k+l)T sample and hold switch 68 providing an odd 
sampled phase Q AQ that is communicated connected to an I 
hard limiter 70. The I hard limiter 70 provides zero and one 
hard limited estimated data A„ representing the odd sampled 
phase 0 A£? from the sampled and held switch 68 sampling at 
odd sample times (2k+l)T. The inphase filtered signal 0, is 
sampled and held by switch 68 at the beginning of the odd 
sample times for one bit period T. That is, the sampling 
40 switch 66 provides the even sample phase 0^ by sampling 
the real PAM filtered signal 0 7 at the beginning of even 
sample times 2kT and holding the sampled value for one bit 
duration T. The 0 G PAM filtered signal from the Imaginary 
PAM filter 64 is communicated to a 2kT sample and hold 
45 switch 72 providing a Q even sampled phase B B/ and to a 
(2k+l)T sample and hold switch 74 providing a Q odd 
sampled phase 0 BQ that is then communicated to a Q hard 
limiter 76. The Q BJ sampled phase output from the sampling 
switch 72 is fed into a Q hard limiter 76 for providing a hard 
50 limited estimated data S„. The estimated data sequences £„ 
and A„ are estimated data of the precoded sequences 
because the absolute phase of the phase of real and imagi- 
nary PAM filler signals 0 7 and Q Q correspond to the precoded 
data sequence a(t) at symbol boundaries. 
55 Sampling switches 66 and 72 operate in synchronism 
clocking at 2kT even symbol boundaries, whereas sampling 
switches 68 and 74 operate in synchronism clocking at 
(2k+l)T odd symbol boundaries. The term k is an incre- 
menting integer index for consecutive symbol boundary 
60 times. The switches 66, 68, 72 and 74 are switched at 
multiples of T, the symbol time period, when the outputs of 
66, 68, 72 and 74 are held through the symbol time T. 
Switches 66 and 72 synchronously respectively sample Q 
and I PAM filtered outputs 0 O and 0 7 at times 2kT, and 
65 switches 68 and 74 sample Q and I PAM filtered outputs d Q 
and 0, at times (2k+l)T. The output of the sampling switches 
66, 68, 72 and 74 provide sampled phase values that are 
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positive and negative and between +1 and -1 indicating a 
phase rotation and therefore indicating binary values of data 
bits of the precoded sequence a(t). 

Positive and negative encoded bits B„ and A„ are 
demodulated, converted, filtered and sampled to provide 
sample phase outputs 0^, 0^, Q BJi and B BQ , from switches 
66, 68, 72 and 74, respectively. The sampled phase outputs 
Q AQ is hard limited by the Q hard limiter 70 to provide the 
real data estimate A„, and the sampled phase outputs 0 BJ is 
hard limited by^the I hard limiter 76 to provide the imaginary 
data estimate B n . The real and imaginary data estimates A„ 
and 6„ provide alternating estimated data bits, that when 
combined, reconstitutes the original data stream d(t). The Q 
and I multipliers 78 and 80 determine the actual phase 
difference between perfect demodulation of the data esti- 
mates fi„ and A„ with the current phase 0^ and Q BQ , 
respectively, to derive phase error signals ^ en and €. odd , 
respectively, that are summed by summer 82 into a total 
phase error e„. 

Switches 66 and 72 are orthogonally switched in synchro- 
nism and operate in combination with hard limiters 70 and 
76, multipliers 78 and 80, and summer 82 to provide a total 
error signal e w . The hard limiter outputs A„ and S n provide 
absolute data signals as sign indicators for determining 
the total error signal e„. The sampled phase outputs 9^ and 
& BQ are values providing a measure of how close the phase 
values are to zero at the time of orthogonal clocking of the 
hard limiters 70, and 76, respectively, generating actual data 
estimates at perfect modulation. The phase output value of 
0^ from switch 66 is at zero during perfect timing when 
hard limiter 76 is orthogonally clocked, and the phase output 
value of d BI is at zero during perfect timing when hard 
limiter 70 is orthogonally clocked. Hence, the hard limiters 
provide alternating data estimates B„ and A„ using orthogo- 
nally sampled filtered signal 0^, and 0^. The sample phase 
outputs 0^ and Q BQ are at zero during perfect timing when 
hard limiters 70 and 76 are orthogonally clocked. When the 
clocking is not perfect, the output of the multipliers 78 and 
80 will provide I and Q phase errors e,.^ and e^, 
respectively, that are summed by the summer 82 providing 
the total error signal e„. 

The switches 66 and 72 synchronously switch at even 
number symbol boundaries. When the switch 66 is clocking 
on even number symbol boundaries, the real inphase phase 
of the I PAM filtered signal 0 7 is near to a zero or small phase 
shift value. The larger the phase shift error, the more out of 
phase the clocking is to the true symbol boundary and hence 
the larger is the real inphase phase shift error 0 M . The 
imaginary quadrature Q PAM filtered signal Q Q is at a 
maximum phase for providing the bit estimate 6„ for 
providing a correct sign of the phase error as sampled by 
switch 66 on the even symbol boundaries. The I multiplier 
78 multiplies the imaginary data estimate 8„ as the sign of 
the real phase error, by the magnitude of the real inphase 
phase error 0^ to obtain a sign and magnitude product 
defining the phase error of the real I PAM filtered signal 0/ 
to provide the inphase phase error e^„. 

The switches 68 and 74 synchronously switch at odd 
number symbol boundaries. When the switch 74 is clocking 
on odd number symbol boundaries, the imaginary quadra- 
ture phase of the Q PAM filtered signal Q Q is near or at a zero 
phase shift value. The larger the phase shift error, the more 
out of phase the clocking is to the true symbol boundary and 
hence the larger is the quadrature phase error Q BQ . The real 
inphase I PAM filtered signal 6 7 is at a maximum phase for 
providing the bit estimate A„ for providing a correct sign of 
the phase error Q BQ as sampled by switch 74 on the odd 



symbol boundaries. The Q multiplier 80 multiplies the sign 
of the real data estimate the imaginary quadrature phase 
error, by the magnitude of the imaginary quadrature phase 
error B BQ to obtain a sign and magnitude product defining 

5 the phase error of the real I PAM filtered signal 0 Q to provide 
the imaginary quadrature phase error e^. 

The switches 66 and 72 clock at even number symbol 
boundaries as the 68 and 74 switches clock at odd number 
symbol boundaries. When the switch 66 is clocked, the real 

10 inphase phase 0 7 is close to zero as a small phase shift. The 
larger and more out of phase is the real inphase phase 6, 
respecting perfect timing, the larger is the imaginary phase 
error when switch 72 is clocked when the imaginary quadra- 
ture filtered signal Q Q is at a maximum phase providing the 

15 bit estimate 6„ providing the correct sign of the real inphase 
phase error. The I multiplier 78 multiplies the imaginary data 
estimate B„ by the real inphase phase 0^ to obtain the sign 
and magnitude of the phase shift of the real inphase phase 
error t even . When the switch 74 is clocked, the imaginary 

20 quadrature phase B Q is close to a zero as a small phase shift. 
The larger and more out of phase is the imaginary quadrature 
phase 0 G respecting perfect timing, the larger is the real 
quadrature phase error when switch 68 is clocked when the 
real inphase phase filtered signal B Q is at a maximum phase 

25 providing the bit estimate A„ providing the correct sign of 
the imaginary quadrature phase error. The Q multiplier 80 
multiplies the real data estimate A„ by the imaginary quadra- 
ture phase Q BQ to obtain the sign and magnitude of the 
imaginary quadrature phase error I multiplier 78 mul- 

30 tiplies the sampled phase shift 0^ with the fi„ hard limited 
estimated data to provide the even phase error e^^, while a 
Q multiplier 80 multiplies the sampled phase shift 0 G with 
the A„ hard limited estimated data to provide the odd phase 
error e^. The even and odd phase errors e tfV(f „ and are 

35 summed by summer 82 to provide the total error signal e„. 
The GMSK signal starts at a zero phase for the first 
precoded symbol. Thereafter, the phase modulation is con- 
tinuous and varies between +/-ji/2. The hard limiters 70 and 
76 provide respective data symbol estimates A„ and B„ at the 

40 symbol time T, and hence functions as first order data 
detectors through a defined relationship between the phase 
errors and the estimated data. 
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cos(e„)+sin(e„); for n=cvcn 
€„**4 cos(€„)+sin(e„); for /i=odd 

A positive one estimate data corresponds to +ji/2 phase 
rotation and a negative one corresponds to -Jt/2 phase 
rotation. At a symbol T, the multipliers 78 and 80, 
respectively, multiply the data estimates B„ and A„ by 
respective orthogonal phase errors 0^ and d BQ . The I 
multiplier 78 multiplies the imaginary quadrature data esti- 
mate 6 n by the real inphase phase error 0^ which is the real 
phase error at xj2 that is typically very small, for determin- 
ing the real inphase even phase error c ewn . Likewise, the Q 
multiplier 30 multiplies the real data estimate A„ by the 
imaginary quadrature phase error Q BQ , which is the real 
inphase phase error at rz/2 that is also typically very small, 
for determining the imaginary quadrature odd phase error 
*-odd- Th e summer 82 adds the phase errors and t even to 
provide the composite total error signal e„. 

The summer 82 receives the imaginary quadrature odd 
phase error from multiplier 80 during odd symbol times 
(2k+l)T when the output from multiplier 78 is zero, assum- 
ing perfect timing, or is receiving the real inphase even 
phase error e even during even symbol times 2kT from the 
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multiplier 78 when the output from the Q multiplier 80 is 
zero, assuming perfect timing. The summer 82 inverts the 
real inphase even phase error e tfWrt for proper summation. 
The real inphase even phase error c even is inverted so that 
when subtracted by summer 82, the real inphase even phase 
error e even is added to the imaginary quadrature odd phase 
error to provide the total error signal e„. The total error 
signal e n has a value of either +1 or -1 for each bit period. 
The loop filter 88 averages the total error value e„ to provide 
a filter error value to the NCO 90 that translates the filter 
error value into a phase error estimate 9 tf . 

The error signal e n is a constant value over a symbol time 
duration T. The absolute value generator 84 desensitizes the 
e„ error signal to noise spikes so that the error estimate is 
constantly nearly one during each symbol tinu period. The 
loop filter 88 digitally filters the error estimate e„ to provide 
a numerically controlled oscillator (NCO) 90 with a stable 
filtered phase error estimate. Trie NCO translates the filtered 
phase error estimate into the estimated phase error 0 tf . The 
estimated phase error 0 tf is added to the carrier phase G c and 
this sum is then used by the phase demodulator 56 to correct 
the carrier phase for accurate closed loop phase demodula- 
tion of the receive signal baseband signal RJi). The com- 
plex receive signal R(t) can be expressed in terms of the real 
and imaginary parts. 

e*-A n +J B n 

The precoder 20 precodes the input data d(t) so that the 
signal magnitude of the filtered signals represent the input 
data d(t) at the odd and even sample times. Hence, the Data 
is present during odd and even samples, and can be sampled 
for improved phase demodulation. The data 6„ and A„ 
sampled by switches 68 and 72 and hard limited by limiters 
70 and 76 and the phases 6 BQ and 6^ sampled by the 
switches 66 and 74 are then used by multipliers 80 and 78 
and summer 82 to generate the phase error signal e rt . The 
error signal e„ is in the form of a phase error. 

e**' = cos(0,) = j srin(0 f ) 

= (A„cos(0 r ) - B n cos(9 r )) +>(*„cos(0 r ) + A a coi@ e )) 
E n = odd = Rete*^*'] = /4„cos(0 f ) - 4? A cos(0 r ) 
e n = even = Im[e> , ^ +fl ' ) ] = S n cos(0 r ) + A„cos(0 r ) 

M„cos(0,) - fl n cos(0,)) + j(B n cos(0 e ) + A n cos(6 e )) 



e n s 2sin(0 r ) 



The present invention is directed to a Laurent filtering 
GMSK carrier tracking loop operating at baseband for 
demodulating a received GMSK signal communicating pre- 
coded data, into a demodulated baseband received signal. 
No squaring of the received signal is required. The invention 
is characterized by separating the demodulating baseband 
receive signal into real and imaginary components both 
subjected to Laurent filtering to provide baseband soft data 
signals of encoded data subject to orthogonal switching and 
hard limiting to generate phase error signals at baseband for 
closed loop carrier phase tracking. The carrier tracking loop 
does provide estimates of the data sequence and could be 
modified to provide data estimates as part of data detection. 
Those skilled in the art can make enhancements, 
improvements, and modifications to the invention, and these 
enhancements, improvements, and modifications may none- 
theless fall within the spirit and scope of the following 
claims. 



What is claimed is: 

1. A method of demodulating a received signal having a 
phase modulated carrier signal of a carrier signal modulated 
by symbols having a symbol time period, the symbols are 

5 generated from preceding an input data stream as a preceded 
input data stream, the method comprising the steps of, 
carrier demodulating the received signal by a carrier 

reference into a carrier demodulated signal, 
phase demodulating the carrier demodulated signal by a 
10 phase amount into a baseband received signal, 

converting the baseband received signal into an imaginary 

component and a real component, 
real filtering the real component into a real phase signal, 
imaginary filtering the imaginary component into an 
15 imaginary phase signal, 

even real sampling at even symbol time boundaries the 
real phase signal into a sampled even real phase signal, 
odd real sampling at odd symbol time boundaries the real 
phase signal into a sampled odd real phase signal, 
20 real hard limiting the sampled odd real phase signal into 
real data estimates, 
even imaginary sampling at even symbol time boundaries 
the imaginary phase signal into a sampled even imagi- 
nary phase signal, 
25 odd imaginary sampling at odd symbol time boundaries 
the imaginary phase signal into a sampled odd imagi- 
nary phase signal, 
imaginary hard limiting the sampled even imaginary 
phase signal into imaginary data estimates, the imagi- 
30 nary and real data estimates are alternating estimates of 
the input data stream, 
even phase multiplying the imaginary data estimates by 
the sampled even real phase signal for generating an 
even phase error, 
35 odd phase multiplying the real data estimates by the 
sampled odd imaginary phase signal for generating an 
odd phase error, 
combining the even and odd phase errors into a phase 
error estimate, and 
40 phase adjusting the phase amount by the phase error 
estimate for phase demodulation of the carrier demodu- 
lated signal during the phase demodulating step. 

2. The method of claim 1 wherein, 
the precoded input data stream is Gaussian filtered to 

provide the received signal having a pulse amplitude 
modulation representation of the Gaussian filtering, 
the real and imaginary filtering steps are matched filtering 
steps that are matched by the pulse amplitude modu- 
lation representation to the Gaussian filtering, and 
the real and imaginary phase signals are filtered signals 
having an absolute phase at a periodic sampling time 
for indicating the data of the input data stream. 

3. The method of claim 1 wherein, 
the precoded input data stream is Gaussian filtered to 

provide the received signal having a pulse amplitude 
modulation representation of the Gaussian filtering, 
the real and imaginary filtering steps are matched filtering 
steps that are matched by the pulse amplitude modu- 
60 lation representation to the Gaussian filtering, 

the real and imaginary phase signals are filtered signals 
having an absolute phase at a periodic sampling time 
for indicating the data of the input data stream, and 
the matched filtering steps apply a principal Laurent 
65 function to the baseband signal so that the real and 
imaginary filtered signal comprises a Laurent compo- 
nent. 
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4. The method of claim 1 wherein, 

the even phase multiplying step determines a phase dif- 
ference between perfect demodulation of the imaginary 
data estimates and an actual demodulation of the imagi- 
nary data estimates, that difference being indicated by 5 
a sign of the sampled real phase signal for determining 
the sign of the even phase error, and 

the odd phase multiplying step determines a phase dif- 
ference between perfect demodulation of the real data 
estimates and an actual demodulation of the real data 10 
estimates, that difference being indicated by a sign of 
the sampled imaginary phase signal for determining the 
sign of the odd phase error. 

5. The method of claim 1 wherein the combining step 
comprises the steps of, 35 

summing the odd phase error and the even phase error into 
a summed phase error, 

absolute error generating an absolute phase error from the 
summed phase error, 20 

dividing the summed phase error by the absolute phase 
error for providing a constant phase error signal having 
a toggling sign, 

loop filtering the constant phase error to provide an error 
estimate, and 25 

generating the phase error estimate from the error esti- 
mate. 

6. The method of claim 1 further comprising the step of, 
receiving a carrier phase estimate, the phase adjusting steps 
adjust the carrier phase estimate by the phase error estimate 30 
for determining the phase amount. 

7. The method of claim 1 further comprising the step of, 
receiving the carrier reference as a carrier reference estimate 
of the carrier signal. 

8. A method of generating an output data stream from a 35 
received signal having a phase modulated carrier signal of a 
carrier signal modulated by precoded symbols having a 
symbol time period, the symbols are generated from pre- 
ceding an input data stream as a precoded input data stream, 

, the method comprising the steps of, 40 
generating a carrier reference estimate from the received 
signal, 

generating a carrier phase estimate from the received 
signal, 45 

carrier demodulating the received signal by the carrier 
reference estimate into a carrier demodulated signal, 

phase demodulating the carrier demodulated signal by a 
phase amount into a baseband received signal, 

converting the baseband received signal into an imaginary 50 
component and a real component, 

real filtering the real component into a real phase signal, 

imaginary filtering the imaginary component into an 
imaginary phase signal, 

even real sampling at even symbol time boundaries the 55 
real phase signal into a sampled even real phase signal, 

odd real sampling at odd symbol time boundaries the real 
phase signal into a sampled odd real phase signal, 

real hard limiting the sampled odd real phase signal into 60 
real data estimates, 

even imaginary sampling at even symbol time boundaries 
the imaginary phase signal into a sampled even imagi- 
nary phase signal, 

odd imaginary sampling at odd symbol time boundaries 65 
the imaginary phase signal into a sampled odd imagi- 
nary phase signal, 
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imaginary hard limiting the sampled even imaginary 
phase signal into imaginary data estimates, the imagi- 
nary and real data estimates are alternating estimates of 
the input data stream, 

even phase multiplying the imaginary data estimates by 
the sampled even real phase signal for generating an 
even phase error, 

odd phase multiplying the real data estimates by the 
sampled odd imaginary phase signal for generating an 
odd phase error, 

combining the even and odd phase errors into a phase 
error estimate, 

phase adjusting the phase amount by the phase error 
estimate for phase demodulation of the carrier demodu- 
lated signal during the phase demodulating step by 
adjusting the carrier phase estimate by the phase error 
estimate for determining the phase amount, and 

data detecting the baseband received signal into the output 
data stream being an estimate of the input data stream. 

9. The method of claim 8 wherein, 

the precoded input data stream is Gaussian filtered to 
provide the received signal having a pulse amplitude 
modulation representation of the Gaussian filtering, 

the real and imaginary filtering steps are matched filtering 
steps that are matched by the pulse amplitude modu- 
lation representation to the Gaussian filtering, and 

the real and imaginary phase signals are filtered signals 
having an absolute phase at a periodic sampling time 
for indicating the data of the input data stream. 

10. The method of claim 8 wherein, 

the precoded input data stream is Gaussian filtered to 
provide the received signal having a pulse amplitude 
modulation representation of the Gaussian filtering, 

the real and imaginary filtering steps are matched filtering 
steps that are matched by the pulse amplitude modu- 
lation representation to the Gaussian filtering, and 

the real and imaginary phase signals are filtered signals 
having an absolute phase at a periodic sampling time 
for indicating the data of the input data stream, and 

the matched filtering steps apply a principal Laurent 
function to the baseband signal so that the real and 
imaginary filtered signal comprises a Laurent compo- 
nent. 

11. The method of claim 8 wherein, 

the even phase multiplying step determines a phase dif- 
ference between perfect demodulation of the imaginary 
data estimates and an actual demodulation of the imagi- 
nary data estimates, that difference being indicated by 
a sign of the sampled real phase signal for determining 
the sign of the even phase error, and 

the odd phase multiplying step determines a phase dif- 
ference between perfect demodulation of the real data 
estimates and an actual demodulation of the real data 
estimates, that difference being indicated by a sign of 
the sampled imaginary phase signal for determining the 
sign of the odd phase error. 

12. The method of claim 8 wherein the combining step 
comprises the steps of, 

summing the odd phase error and the even phase error into 

a summed phase error, 
absolute error generating an absolute phase error from the 

summed phase error, 
dividing the summed phase error by the absolute phase 

error for providing a constant phase error signal having 

a toggling sign, 
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loop filtering the constant phase error to provide an error 
estimate, and 

generating the phase error estimate from the error esti- 
mate. 

13. The method of claim 8 further comprising the steps of, 
generating a timing estimate from the received signal, and 
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generating a timing signal from the timing estimate and 
the baseband received signal, the timing estimate is a 
clocking signal for sampling the baseband received 
signal during the detecting step. 
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